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ABSTRACT 



An experimental investigation was conducted of the effects 
of fluid loading on axially symmetric flexural waves in a 
circular aluminum plate 25 in. in diameter and 5/16 in. thick. 
Measurements were made of the flexural wavelength and natural 
frequencies with the plate in air, and with water loading on 
one side, over a frequency range from 20 kHz to 80 kHz. The 
critical (coincidence) frequency, that frequency at which the 
speed of the flexural wave equals the speed of sound in water, 
occurs at approximately 45 kHz. The measured flexural-wave 
speeds with air loading are in good agreement with theory. 

With water loading, the wave speed decreases, the amount of 
decrease increasing to approximately 14 percent below the in- 
air value just below the critical frequency. At the critical 
frequency, the wave speed appears to jump to values nearly 
equal to the unloaded values. Radiation patterns were ob- 
tained and when edge effects are eliminated, are in qualita- 
tive agreement with theoretical predictions for an infinite 
fluid- loaded plate. 
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I. INTRODUCTION 



The purpose of this study was to provide some basic infor- 
mation which would support the problem of designing and devel- 
oping a f lexural-wave transducer suitable for the high frequency 
(75 kHz) acoustic ranging system used at the underwater weapons 
test ranges of the Naval Torpedo Station, Keyport, Washington. 
Previous work on the subject [Ref. 1] indicated the possibility 
of using a flexural -wave transducer as a broad-beam, broad- 
bandwidth sound source which would be readily adaptable to 
flush mounting, a necessary requirement for high speed torpedos. 
The possibility of using the torpedo hull itself as the radiat- 
ing surface is especially attractive for low-drag, laminar flow 
shapes . 

In the report of his experimental work on a flexural disk 
transducer, Sevdik [Ref. 1] discusses discrepancies between 
theoretical normal mode frequencies and measured values which 
were obtained for his transducers. He also discussed damping 
of the plate natural frequencies by fluid loading and incon- 
sistencies in the change of phase speed with frequency for 
transducers of differing thickness. 

Sevdik based his calculations on classical plate theory as 
described in Morse and Ingard [Ref. 2] and Malecki [Ref. 3] . 
Mindlin [Ref. 4] includes the effects of shear and rotatory 
inertia in a three-dimensional plate theory and states that 
classical plate theory begins to differ markedly from 
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three-dimensional theory at frequencies above that at which 
the flexural wavelength is about ten times the plate thickness. 
Mindlin's theory gives flexural-wave speed values for a given 
ratio of plate thickness to flexural wavelength which are 
lower than those from classical plate theory. 

Alper and Magrab [Ref. 5] extend the Mindlin-Timoshenko 
plate theory to fluid-loaded plates and present calculated 
results showing decreasing plate natural frequencies resulting 
from fluid loading. 

Feit [Ref. 6] considers the case of an infinite plate which 
is fluid loaded on one side and discusses significant changes 
in the theoretical radiation patterns that result at frequen- 
cies above and below the critical frequency, that frequency 
at which the flexural- wave speed in the plate equals the speed 
of sound in the fluid. Feit's work is also based on the 
Mindlin-Timoshenko theory. 

The objective has been to determine experimentally the 
effects of fluid loading on the dispersion curves for flexural 
waves in a plate. The methods used included observation of 
the natural frequencies in a circular, free-edge plate in air 
and in water, and in .some cases, measurements of the wave- 
lengths of the flexural waves. The frequency range bracketed 
that for the critical wave speed. 

A secondary object was to measure the radiation patterns 
for this circular plate. Some preliminary results are 
reported. 
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II. THEORETICAL CONSIDERATIONS 



A. CLASSICAL PLATE THEORY 

The classical differential equation of motion for the 
transverse displacement W of a plate is given by [Ref. 7, 



1] 



v * w+ T 



d 2 w 

at 2 



- o 



Flexural wave speed is given by [Ref. 3, p. 508] 

./oiV* 



=(£ 



■J 



( 2 . 1 ) 



( 2 . 2 ) 



Classical p late- vibration theory does not include the effects 
of rotary inertia, which were shown to reduce the natural 
frequencies of beams by Lord Rayleigh [Ref. 8] , nor the effects 
of shear deformation, which were shown to be equally important 
by Timoshenko [Ref. 9]. The classical theory is valid for the 
case where the plate is thin compared to the flexural wave- 
length. 



B. MINDL IN -TIMOSHENKO PLATE THEORY 

Mindlin [Ref. 4] extended Timoshenko's theory for beams 
by including the effects of both shear and rotatory inertia 
in a three-dimensional plate theory. He has shown that the 
shear correction accounts almost entirely for the difference 
between classical plate theory and three-dimensional plate 
theory over the whole wavelength spectrum, and derived rela- 
tionships for flexural wave speed as a function of the ratio 
plate thickness/flexural wavelength. For the case where only 
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the shear correction is considered [Ref. 4, p. 37] 
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Substituting for c ^ in equation 2.3 and rearranging 
gives 
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Equation 2.5 reduces to 
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Starting with the classical plate theory equation for flexural 
wave speed (2.2), substituting 9 h for m, substituting for D, 



and squaring gives 
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Substituting oj = 

C =A(tf 
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and rearranging gives 



(2.9) 



( 2v 10) 



Equations (.2.8) and (2.10) relate flexural wave speed to 
for Mindlin’s shear-corrected theory and classical plate theory, 
respectively. The constants A and B are determined by the 
properties of the plate material. For aluminum [Ref. 10]: 
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p = 2700 kg/m 3 
E = 7.1 x 10 10 Nt/m 2 
q = Q.33 
K 2 = 0.8686 
A = 9.71 x 1 0 10 (m/ s) 2 
B = 11.3 

Using these values, flexural-wave speed is plotted as a 
function of h/A in Fig. 1. It can be seen that the shear- 
corrected value for flexural- wave speed begins to differ 
significantly from that given by classical plate theory for 
values of j > 0.10. For Sevdik's 3 mm thick transducer at 
at 75 kHz, j = 0.15. 

C. EFFECTS OF FLUID LOADING 

Feit [Ref. 6] considers the case of an infinite, elastic 
plate, fluid-loaded on one side and excited by a time-harmonic 
point force. He derives an expression for the acoustic power 
radiated by this plate into an infinite half-space filled with 
an acoustic fluid. Because the velocity of sound in water is 
relatively high compared to the velocity of sound in air, the 
frequency at which the flexural wave speed in the plate exceeds 
the speed of sound in the fluid, i.e., the critical frequency, 
is relatively high. For this reason, Feit uses the Mindlin- 
Timoshenko plate theory in his work. Feit's equation for the 
farfield pressure is [Ref. 6, p. 1492] 
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PLATE THICKNESS/FLEXURAL WAVELENGTH 
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and his expression for the critical frequency is 




Feit’s theory predicts significant changes in the radiation 

patterns from the infinite plate at frequencies above and 

below the critical frequency. Quoting Feit [Ref. 6, p. 1493], 

Below coincidence, the propagating flexural wave does 
not radiate energy away from the plate. Therefore, the 
radiated sound pressure can be associated with the non- 
propagating nearfield deflection of the plate, which 
acts like a small piston and therefore is relatively non- 
directive . 

Above coincidence, the situation is quite different. 

In this range, the phase speed of the propagating flexural 
wave exceeds that of the sound wave, and therefore the 
propagating wave can itself radiate energy away from the 
plate in a definite direction. 

Alper and Magrab [Ref. 5] extend the Mindlin-Timoshenko 
plate theory to the case of clamped-edge circular plates with 
fluid loading on one side. They present results based on 
computations using their theory for the case of a thin 
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(thickness/radius = 0.04) and a thick (thickness/radius = 0.2) 
steel plate in sea water. Their results indicate a signifi- 
cant reduction in the plate natural frequencies when fluid 
loaded which is more pronounced for the thin plate. They also 
conclude that for the same driving force the thin plate will 
radiate considerably more power into the water and that the 
thin plate produces a more directional radiation pattern for 
a given driving frequency. 
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III. EXPERIMENTAL PROCEDURE 



A. FLEXURAL WAVE SPEED MEASUREMENTS 

Several experimental techniques were used to measure 

flexural wave speed as a function of frequency. 

Initial measurements were made with a 1/4 in. thick 

rectangular steel plate, 64.5 cm by 76.5 cm. The plate was 

driven by a longitudinal vibrator composed of four 1/4 in. 

thick, 3/4 in. diameter piezoceramic elements. Brass masses 

1/4 in. thick and 3/4 in. in diameter were added at each end 

of the ceramic elements and the unit was clamped together 

with a screw threaded into the bottom disk. The unmounted 

driver was resonant at 34.01 kHz. The driver was mounted 

near the geometric center of the plate with epoxy cement. 

With the driver mounted on the plate, over 105 resonances 

were detected between 10 kHz and 40 kHz using measurements 

of the electrical admittance of the driver transducer. 

Two different experimental techniques were employed in 

measuring flexural wave speed. The first procedure was to 

use a small accelerometer to measure the phase and amplitude 

of the wave in the plate along a radial line from the driver 

transducer while the plate was excited continuously at one 

frequency. The equipment used in this procedure is listed 

below and connected as shown in Fig. 2. 

Hewlett Packard 3590A/3594A Wave Analyzer/Oscillator 

Dranetz Model 305 Digital Phase Meter 

Tektronix 547 Dual Trace Oscilloscope 

Hewlett Packard 465A Voltage Amplifier 

Hewlett Packard 467A Power Amplifier 

Data Precision 5740 Frequency Counter 
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Figure 2. Phase/Amplitude Measurement Apparatus 
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Using the driving voltage as the reference for the phase- 
meter and the voltage generated by the accelerometer to 
measure the phase of the motion of the plate, measurements 
were made at one-half cm intervals over a distance of 15 cm. 
From the measured values of total phase change over the known 
distance, flexural wavelength could be determined. 

Two difficulties were encountered using this technique. 
First, acquiring data was a slow and time-consuming process. 

A thin layer of grease was smeared along the line over which 
the measurements were made to get good contact between the 
accelerometer and the plate. It was found that the values 
of phase and amplitude obtained would vary depending upon how 
the accelerometer was placed on the plate. To compensate for 
this variation and positioning error, a minimum of three sets 
of measurements were made at each frequency. This required 
a minimum of 90 measurements to obtain one value of flexural 
wavelength versus frequency. Secondly, it was found that 
because of the very complicated modal structure in the plate 
due to asymmetry, the phase data did not give reliable wave- 
length values. Due to these difficulties, this technique was 
abandoned. 

The second technique used was to excite the plate with a 
short (four cycle) pulse. By measuring the travel time of a 
specific cycle of the pulse over a known distance, flexural 
wave speed could be determined. This method also eliminated 
the complex standing wave patterns resulting from continuous 
excitation. The test equipment used for these measurements 
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is listed below and depicted schematically in Fig. 3. 

Hewlett Packard 3590A/3594A Wave Analyzer/Oscillator 

Data Precision 5740 Frequency Counter 

General Radio 1398A Tone Burst Generator 

Tektronix 547 Dual Trace Oscilloscope 

Hewlett Packard 465A Voltage Amplifier 

Hewlett Packard 467A Power Amplifier 

40 dB Gate/Harrison 6205B DC Power Supply 

The driving voltage was used as a reference and wave 
travel time measurements were made using the delayed sweep 
feature of the oscilloscope. Six measurements were made 
at each frequency and averaged to compensate for error in 
positioning the accelerometer. These data are plotted in 
Fig. 4 in comparison with a curve obtained from Mindlin 
(Eqn. 2.8). Error bars, the average deviation above and 
below the mean of the six values, are plotted to indicate 
typical variability of the data. 

The next procedure was to measure the flexural- wave speed 
with one surface of the plate in contact with water. A three- 
inch high plywood dam, two feet on a side, was secured and 
sealed to the plate which permitted making measurements on 
the top of the plate with the accelerometer while the plate 
essentially floated horizontally in the anechoic tank. The 
pulse method was used to measure wave speed and six measure- 
ments were made and averaged, as was done with the in- air 
data. These data are also plotted with the in-air data in 
Fig. 4. As can be seen, the in-water data do not differ sig- 
nificantly from the in-air data. The reason for this result 
is not certain since it was expected that the flexural- wave 
speed values would be lower for a given frequency in water. 
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Figure 3. Pulse Measurement Apparatus 
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Figure 4. Steel Plate Flexural Wave Speed 



Work with the steel plate was discontinued at this point due 
to the complexity of the modal structure and the difficulty 
in handling the large weight. 

Measurementswere continued on a round aluminum plate, 

25 in. in diameter and 5/16 in. thick. It was driven at the 
center by the same stack of ceramic elements and a three-inch 
high sheet aluminum dam was secured and sealed to the perimeter 
of the plate to permit making measurements with one side of 
the plate in contact with water. The plate with dam was 
buoyant and floated approximately one inch deep in the anechoic 
tank. Using the pulse method, six sets of data were obtained 
in air and in water. These data were even more irregular and 
scattered than the data obtained with the steel plate. How- 
ever, it was observed that dispersion effects in the propaga- 
tion of the wave pulse contributed significantly to the 
problem of measurement of the phase speed. These effects may 
also have been a contributor to the problems with measurements 
on the steel plate described above. Because of the difficulties 
encountered using the pulse method to obtain data, it was 
discontinued. 

B. PLATE RESONANCE FREQUENCY: MEASUREMENTS 

Another approach to the observation of the effects of fluid 
loading on the flexural waves in the plate was used. This was 
accomplished by measuring the frequencies of the f lexural-wave 
normal modes for the free-edge circular plate. The electrical 
admittance of the driving transducer provides one method for 
observing the plate resonances. The equipment listed below 
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was connected as shown in Fig. 5 for this procedure: 

Hewlett Packard 3590A/3594A Wave Analyzer/Oscillator 

Hewlett Packard 7035B X-Y Plotter 

Hewlett Packard 465A Voltage Amplifier 

Hewlett Packard 467A Power Amplifier 

Dranetz 100C Complex Impedance/Admittance Meter 

Data Precision 5740 Frequency Counter 

As can be seen from the equipment diagram, the wave analyzer 
was used to drive the X axis of the plotter. Plots of con- 
ductance versus frequency were made rather than plots of 
conductance and susceptance because they presented a more 
graphic picture of the normal mode frequencies. Fig. 6 is 
representative of the plots obtained using this technique. 

In an attempt to identify each mode, theoretical values were 
calculated using eigenvalues obtained from Leissa [Ref. 7, 
p. 11]. Because Leissa's values are computed from classical 
plate theory and because of the effect of the mass loading of 
the driver on the plate, the theoretical frequencies could 
not be directly correlated with the experimental values. 
Therefore, an experimental procedure was needed to determine 
which mode was excited at a given frequency. This was carried 
out by mapping the modal structure in the plate using an 
accelerometer. It was found in making measurements at those 
frequencies where peaks occurred in the conductance vs. fre- 
quency plots that only modes with nodal circles were detected. 
The amplitude of the signal generated by the accelerometer 
remained constant when the accelerometer was moved in a circle 
about the center of the plate. It was also found that the 
nodal circles (amplitude minimums) were quite regularly spaced 
for the higher order modes (n > 8) . By measuring the spacing 
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Figure 5. Nodal Circle Measurement Apparatus 
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FREQUENCY (kHz) 

Figure 6. Conductance/Frequency Plot, Aluminum Plate in Air 



between nodal circles at the plate natural frequencies 
detected with the Dranetz meter, flexural wavelength could 
be measured quite accurately as a function of frequency. The 
assumption that the spacing between two nodal circles is equal 
to one-half flexural wavelength is justified in Appendix A. 

The initial nodal circle measurements were made with the 
sheet aluminum dam removed to more closely approach a free 
edge boundary condition. To permit making measurements in 
water, a dam was made from 0.011 in. thick natural rubber 
sheet. The rubber sheet was stretched across the bottom of 
the plate and secured to the sheet metal ring used previously, 
with about one inch of free rubber between the plate and the 
supporting ring. The plate was suspended using heavy nylon 
cord fastened at eight regularly spaced points around the top 
of the aluminum ring. To ensure the plate was loaded by the 
water, the thin space between the bottom of the plate and the 
rubber sheet was filled with transducer oil having a specific 
acoustic impedance equal to that of water and all air bubbles 
carefully removed. This arrangement is depicted in Fig. 7. 

The stack of four ceramic elements used in the previous 
experiments was used to excite the plate. Fig. 8 is the con- 
ductance vs. frequency plot for the plate suspended in the 
rubber dam in air. The flexural wavelength data obtained 
using the nodal-circle method is summarized in Table I and 
plotted in Fig. 9, Figures 10 and 11 present the conductance 
vs. frequency plots obtained with the plate in water. The 
in-water data are summarized in Table II and are plotted with 
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Figure 7. Plate Suspension Apparatus (not to scale) 
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Figure 8. Conductance/Frequency Plot, Plate in Air in Rubber Dam. 



TABLE I 



NODAL CIRCLE DATA IN AIR 



Frequency 


Radial 

Distance 


Number of 
Wavelengths 


Flexural 

Wavelength 


Flexural 
Wave Speed 


7764 kHz 


18.8 cm 


2 


9.40 cm 


730 m/s 


10048 


25.5 


3 


8.50 


854 


12586 


21.8 


3 


7.30 


919 


15366 


27.5 


4 


6.90 


1060 


18274 


24.4 


4 


6.10 


1115 


21229 


22.4 


4 


5.60 


1189 


24085 


26.0 


5 


5. 20 


1252 


29695 


28.2 


6 


4.70 


1396 


33216 


31.5 


7 


4. 50 


1495 


36794 


24.8 


6 


4.13 


1520 


39983 


23.2 


6 


3. 87 


1547 


46749 


25.2 


7 


3.60 


1683 


51236 


23.9 


7 


3.41 


1749 
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Figure 9. Aluminum Plate Flexural Wave Speed 
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FREQUENCY (kHz) 

Figure 10. Conductance/Frequency Plot, Plate Suspended in Water. 
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FREQUENCY (kHz) 

Figure 11. Conductance/Frequency Plot, Plate Suspended in Water 



TABLE II 



NODAL CIRCLE DATA IN WATER 



Frequency 


Radial 

Distance 


Number of 
Wavelengths 


Flexural 

Wavelength 


Flexural 
Wave Speed 


12164 Hz 


27.2 cm 


4 


6.80 cm 


827 m/s 


14624 


27.7 


4.5 


6.16 


900 


17055 


28.4 


5 


5.68 


969 


19770 


28.4 


5.5 


5.16 


1021 


22293 


29.6 


6 


4.93 


1100 


24978 


27.6 


6 


4.60 


1149 


27770 


26. 2 


6 


4.37 


1213 


30604 


24.4 


6 


4.07 


1245 


33718 


23.3 


6 


3.88 


1309 


37155 


22.0 


6 


3.67 


1362 


39633 


20.9 


6 


3.48 


1381 


42827 


19.4 


6 


3.23 


1385 


46002 


18.6 


6 


3.10 


1426 


SET 2 


12178 


20.4 


3 


6.80 


828 


14614 


25. 5 


4 


6.38 


932 


17194 


23.2 


4 


5.80 


997 


19775 


26.8 


5 


5.36 


1060 


22329 


25.0 


5 


5. 00 


1116 


25002 


23.5 


5 


4.70 


1175 


27760 


26. 2 


6 


4.37 


1213 


30636 


24.4 


6 


4.07 


1247 


33673 


23.6 


6 


3.93 


1323 


36745 


22.0 


6 


3.67 


1349 


39516 


2Q . 6 


6 


3.43 


1355 


42425 


19.8 


6 


3.30 


1400 
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the in-air data in Fig, 9. As can be seen, the in-air data 
are in good agreement with the Mindlin theory and the flexural 
wave speed with water ’ loading is significantly lower than it 
is in air. The highest wave speed values obtained in water 
were just above 1400 m/s , which is below the sound speed in 
water. Since it was desired to investigate the flexural-wave 
speed in water at frequencies where it exceeded the sound 
speed in water, a new driver was fabricated from two 3/8 in. 
diameter, 1/4 in. thick piezoceramic elements. With this 
driver mounted on the plate, in-air nodal-circle data were 
obtained to above 80 kHz. These data are summarized in Table 
III and plotted with the data obtained previously in Fig. 13. 
Fig. 12 is the conductance vs. frequency plot with the plate 
suspended in air using the small driver, and Fig. 14 is the 
corresponding plot with the plate in water. It exhibits none 
of the structure of the previous plots and in-water data could 
not be obtained using the nodal-circle technique. 

Because the plate natural frequencies could not be detected 
using the small driver with the plate in water, a second tech- 
nique was used to extend the in-water measurements to 80 kHz. 
Using an oscilloscope with vertical and horizontal sweeps, 
Lissajous patterns were used to measure the rate of change 
of the phase difference between the driving voltage and the 
voltage generated by the accelerometer as a function of the 
distance from the driving transducer. These data are summar- 
ized in Table IV. As can be seen, the data obtained by this 
technique differ from the data obtained using the nodal 
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TABLE III 



IN-AIR DATA/NODAL CIRCLE TECHNIQUE WITH SMALL DRIVER 



Frequency 


Radial 

Distance 


Number of 
Wavelengths 


Flexural 

Wavelength 


Flexural 
Wave Speed 


45643 Hz 


21.8 cm 


6 


3.63 cm 


1658 m/s 


49832 


20.8 


6 


3.47 


1728 


54089 


19.7 


6 


3.28 


1776 


58466 


19.4 


6 


3.23 


1890 


62870 


18.0 


6 


3.00 


1886 


67381 


17.1 


6 


2.85 


1920 


71991 


16.9 


6 


2.82 


2027 


76560 


15. 7 


6 


2.62 


2003 


80976 


15.3 


6 


2. 55 


2060 



IN- WATER DATA/ SMALL 

c Radial 

Frequency Distance 


TABLE IV 

DRIVER/ LI SSAJOUS 

Number of 
Wavelengths 


PATTERN TECHNIQUE 

Flexural Flexural 

Wavelength Wave Speed 


40000 Hz 


21.0 cm 


6 


3. 50 cm 


1400 m/ s 


45000 


18.2 


5 


3.64 


1638 


45760 


18.6 


6 


3.10 


1418 


50000 


13.8 


4 


3.45 


1725 


55000 


16.0 


5 


3. 20 


1760 


60000 


15.5 


5 


3.10 


1860 


65000 


18.2 


6 


3. 03 


1971 


74887 


13.2 


5 


2.64 


1977 
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FREQUENCY (kHz) 

Figure 12. Conductance/Frequency Plot, Small Driver, Plate in Air 
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Figure 13. Aluminum Plate Flexural Wave Speed 
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Figure 14. Conductance/Frequency Plot, Small Driver, Plate in Water 



circle technique. The reasons for this observed difference 
in flexural -wave speed values may involve the differing 
measurement techniques. In order to compare the techniques, 
six measurements were made in air using each method between 
28 kHz and 50 kHz over a distance of five wavelengths. The 
values obtained by the two methods agree within five percent 
whereas the observed change in flexural wave speed was 15 per- 
cent. The effect of radiation impedance is clearly changing 
since the normal modes are highly damped above the critical 
frequency, so the observed changes in flexural wave speed may 
be real. 

C. RADIATION PATTERN MEASUREMENTS 

Alper and Mabrab [Ref. 5] present theoretical radiation 
patterns for two cases of steel plates with clamped edges. 

The radiation patterns are computed for frequencies only up 
to the fifth in-fluid resonant frequency which corresponds 
to approximately four kHz for the circular aluminum plate 
used in this work. 

Feit [Ref. 6] presents theoretical results for an infinite 
plate which is fluid loaded on one side. Based on the Mindlin- 
Timoshenko theory, he predicts significant differences in the 
radiation patterns for the infinite plate at frequencies above 
and below the critical frequency. Below the critical frequency 
the theoretical beam patterns are nearly omnidirectional, 
whereas above the critical frequency theory predicts the 
radiation pattern to be composed of two major side lobes. 

For a frequency equal to 1.5 times the critical frequency, 



41 



Feit's theory predicts these lobes to be approximately 60 
degrees off the plate axis. For a steel plate, Feit’s theory 
predicts that the critical frequency based on the Mindlin- 
Timoshenko theory is approximately 23 percent higher than 
predicted by classical plate theory. 

A series of radiation patterns was made with the circular 
aluminum plate in order to compare theory with experiment. 

The radiation patterns were made in one of the anechoic water 
tanks at the Naval Postgraduate School. 

In order to limit acoustic loading to one side of the 
plate, the back of the plate was covered with a sheet of 1/4 
in. thick neoprene foam, the same type of material used for a 
skin diver's wet suit. The driver transducer was enclosed in 
a plastic cup sealed to the neoprene, which provided electrical 
and mechanical isolation from the water. The neoprene layer 
was glued on at the edge of the plate. 

The below listed equipment was used to make the radiation 

pattern measurements and is depicted schematically in Fig. 15. 

Hewlett Packard 7035B X-Y Plotter 
Tektronix RM 561A Dual Trace Oscilloscope 
Hewlett Packard 3594A Sweeping Oscillator 
General Radio 1396A Tone Burst Generator 
40 dB Gate 

Harrison 6205B DC Power Supply 
Hewlett Packard 465A Power Amplifier (2) 

Krohn-Hite DCA 50R 50-Watt Power Supply 
Hewlett Packard 465A Voltage Amplifier 
Arenberg Ultrasonic Laboratory Attenuator 
Envelope Detector 
PAR Model 160 Boxcar Integrator 

Hewlett Packard/Moseley 7561A Logarithmic Converter 
Data Precision 5740 Frequency Counter 
Hewlett Packard 400E AC Voltmeter 
LC-10 Hydrophone 
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Figure 15. Radiation Pattern Measurement Apparatus 
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The mechanical drive used to rotate the plate includes a 
potentiometer on the gear train which produces a DC voltage 
proportional to the angular position of the plate. This 
voltage was applied to the X axis of the X-Y plotter. The 
calibrated LC-10 hydrophone used as a receiver and the plate 
were separated by a distance of 5 meters with the hydrophone 
and the center of the plate at a depth of 1.92 meters. 

The output of the hydrophone was amplified, rectified by 
the envelope detector and then fed to the boxcar integrator. 
By using a pulsed signal and gating the integrator to inte- 
grate only over the direct path signal, interference effects 
from the air surface reflected signal are eliminated. The 
integrator produces a DC voltage which is proportional to the 
average amplitude of the received pulse envelope. Using the 
logarithmic converter permits the plotting of a large dynamic 
range of received signal. 

The attenuator shown in Fig. 15 was used to calibrate the 
Y axis of the radiation pattern plots. This was necessary 
because of the non-linearity of the envelope detector. Each 
radiation pattern was calibrated by varying the drive voltage 
to obtain a full scale reading on the integrator (10 V. DC 
output from the integrator) and then incrementing the attenu- 
ation in 3 dB intervals. The drive voltage was then reset to 
80 V, RMS prior to generating the radiation pattern. The 
radiation patterns are shown in Figs. 16 through 31, and 
Table V is a summary of the equipment settings for each 
pattern. 
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Figure ]7. Radiation Pattern - 35 kHz/32 Cycle Pulse 
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Figure 19. Radiation Pattern - 45 kHz/32 Cycle Pulse 
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Figure 20. Radiation Pattern - 50 kHz/32 Cycle Pulse 
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Figure 21. Radiation Pattern - 55 kHz/32 Cycle Pulse 
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Figure 22. Radiation Pattern - 60 kHz/32 Cycle Pulse 
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Figure 23. Radiation Pattern - 65 kHz/32 Cycle Pulse 
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Figure 24. Radiation Pattern - 70 kHz/32 Cycle Pulse 
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Figure 25. Radiation Pattern - 30 kHz/4 Cycle Pulse 
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Figure 26. Radiation Pattern - 35 kHz/4 Cycle Pulse 



-12 





56 



DEGREES OFF MAIN AXIS 

Figure 27. Radiation Pattern - 40 kHz/4 Cycle Pulse 
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Figure 28. Radiation Pattern - 45 kHz/4 Cycle Pulse 
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Figure 29. Radiation Pattern - 50 kHz/4 Cycle Pulse 
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Figure 31. Radiation Pattern - 60 kHz/4 Cycle Pulse 



BEAM PATTERN MEASUREMENT EQUIPMENT PARAMETERS 
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The first set of radiation patterns was made using a 
32-cycle pulse. This was the longest pulse that could be 
used without having the water-air surface reflected signal 
interfere with the direct path signal at 30 kHz. The boxcar 
integrator gatewidth was maintained at 500 U sec. The only 
discernable feature from this first set of radiation patterns 
is the highly variable amplitude. By making nearfield meas- 
urements at a distance of one centimeter from the plate, it 
was found that the plate radiated more strongly from the 
center region of the plate and from the edge region of the 
plate. This effect, combined with the large radius (31.75 cm) 
of the plate, resulted in the large amplitude variations as 
the plate rotated. To eliminate the interference effects of 
the waves arriving from the edge region of the plate, a short 
(4 cycle) pulse was used. This permitted separating the waves 
arriving from the central portion and edge region of the plate 
by properly gating the boxcar integrator. By reducing the 
edge effects, a closer approximation was obtained to the 
infinite plate used in Feit's theory. As can be seen, most 
of the interference effects are eliminated and the radiation 
patterns are much smoother. Starting at 45 kHz (Fig. 28), 
the radiation patterns change from a relatively flat, broad, 
single lobe to two side lobes approximately 60 degrees off 
the plate axis at 60 kHz. This is in qualitative agreement 
with Feit's calculations for the infinite fluid-loaded plate. 
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IV. DISCUSSION OF RESULTS 



Because all of the experimental work was done with a 
single plate, it was impossible to arrive at general con- 
clusions. For the data obtained in air, good agreement was 
obtained with the Mindlin-Timoshenko plate theory for the 
speed of the flexural wave as a function of frequency. In 
comparing the in-water data with the in-air data, the flexural- 
wave speed with water loading is approximately 11 percent lower 
than with air loading at 15 kHz and approximately 14 percent 
lower at 45 kHz. Plate natural frequencies were 18 percent 
lower for both the tenth and eighteenth axisymmetric modes, 
occurring in air at 14931 Hz and 40949 Hz, respectively. 

These results agree qualitatively with theory. Additionally, 
the changes which occur in the radiation patterns above 45 kHz, 
which is approximately the critical frequency for the plate, 
appear to agree qualitatively with Feit's theory for the 
infinite fluid-loaded plate. 

The most significant effect noted, which is not predicted 
by any of the theory reviewed by the author, is the behavior 
of the flexural-wave speed above the critical frequency. Al- 
though different measurement techniques were used to obtain 
the data, there is. an indication of a significant (15 percent) 
increase in flexural wave speed occurring in the vicinity of 
45 kHz. It is in this frequency range that the experimental 
data for flexural- wave speed approaches the speed of sound 
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in water. Although the possibility exists that different 
experimental techniques caused this result, the two methods 
of obtaining data agreed within five percent at frequencies 
below the critical frequency and repeated observations indi- 
cate that the result is real. There is clearly a change in 
the radiation resistance as frequency increases, as evidenced 
by the increased damping of the plate natural frequencies 
with increasing frequency. At frequencies below the critical 
frequency, the fluid adds inertia which causes a lower 
flexural wave speed. If this effect is real, it indicates 
that the radiation reactance (inertial) becomes much smaller 
above the critical frequency. 
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V. RECOMMENDATIONS 



The results of this work and those of Sevdik are a promis- 
ing indication of the feasibility of using a flexural disk 
transducer as a high-frequency, broad-beam sound source. 

Other advantages of a flexural disk transducer are its sim- 
plicity of construction, ruggedness, and adaptability for 
flush mounting or as part of the hull structure of the vehicle 
in which it is used. 

An in-depth theoretical study leading to a working model 
transducer, following the theory contained in the Alper-Magrab 
paper, is recommended. The effects of shear and rotatory in- 
ertia are significant at the frequencies in question and must 
be considered. This work and the theory presented by Feit 
indicate that significant changes in radiation patterns do 
occur at frequencies above the critical frequency. Further 
investigation into these changes and the possible changes in 
flexural wave speed observed in this work are recommended. 



65 



APPENDIX A 



COMPARISON OF NODAL CIRCLE SPACING 
WITH FLEXURAL WAVELENGTH 



From Leissa [Ref. 7, p. 11], substituting ft 2 for Leissa's A 2 , 






(A - 1) 



From classical plate theory [Ref. 4, p. 44], 




Rearranging (A- 2) 



Substituting for D in (A-l) using (A.- 5) 



(A- 2) 



(A- 3) 




27TCVQ? 

CA 



(A- 4) 



Substituting w = 2irf and c = fA into (A- 4) 



n z _ i t-TT^a 1 

A* 



Rearranging (A- 5) 

A 27 r 



a 




(A- 5) 



(A- 6) 



Where ft 2 is obtained from Leissa [Ref. 7, Table 2.5]. 

From Leissa [Ref. 7, Table 2.7], the spacing between nodal 
circles for a free-edge, circular plate with Poisson's ratio 
equal to 0.33 equals 0.200 r/a, for the 5,0 mode, where r is 
the distance from the center of the plate to the nodal circle 
and a is the plate radius. 
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Assuming the spacing between nodal circles is one-half 
flexural wavelength, for the 5,0 mode, X = 0.400a. Using 
ft 2 = 245.9 (Table 2.5 for the 5,0 mode) to substitute into 
equation A- 6, 

X =■ -Il£ — = 0.40.1a 
i/245.9 

For the 5,0 mode, classical plate theory gives a value of 
X = 0.401a when solving from calculated eigenvalues and 
0.400 when solving using nodal circle spacing. 
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